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Abstract 

In quantum field theories divergences generally turn up in loop calculations. 
Renormalization is a part of the theory which can be performed only with a proper 
regularization. In low energy effective theories there is a natural cutoff with well 
defined physical meaning, but the naive cutoff regularization is unsatisfactory. 
A Lorentz and gauge symmetry preserving regularization method is discussed in 
four dimension based on momentum cutoff. First we give an overview of various 
regularization methods then the new regularization is introduced. We use the 
conditions of gauge invariance or equivalently the freedom of shift of the loop 
momentum to define the evaluation of the terms carrying even number of Lorentz 
indices, e.g. proportional to The remaining scalar integrals are calculated 

with a four dimensional momentum cutoff. The hnite terms (independent of 
the cutoff) are free of ambiguities coming from subtractions in non-trivial cases. 
Finite parts of the result are equal with the results of dimensional regularization. 
The proposed method can be applied to various physical processes where the use 
of dimensional regularization is subtle or a physical cutoff is present. As a famous 
example it is shown that the triangle anomaly can be calculated unambiguously 
with this new improved cutoff. The anticommutator of 7^ and 7^ multiplied by 
five 7 matrices is proportional to terms that do not vanish under a divergent 
loop-momentum integral, but cancel otherwise. 


1 Introduction 

Several regularization methods are known and used in quantum field theory: three and 
four dimensional momentum cutoff, Pauli-Villars type, dimensional regularization, lat¬ 
tice regularization, Schwinger’s proper time method and others directly linked to renor¬ 
malization like differential renormalization. Dimensional regularization (DREG) [1] is 
the most popular and most appreciated as it respects the gauge and Lorentz symme¬ 
tries. However DREG is not useful in all cases, for example it is not directly applicable 
to supersymmetric gauge theories as it modifies the number of bosons and fermions dif¬ 
ferently. DREG gets rid of (does not identify) naive quadratic divergences, which may 
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be important in low energy effective theories or in the Wilson’s renormalization group 
method. Another shortcoming is that together with (modihed) minimal subtraction 
DREG is a “mass independent” scheme, particle thresholds and decoupling are put in 
the theory by hand [2]. The choice of the ultraviolet regulator always depends on the 
problem. 

In low energy effective held theories there is an explicit cutoff, with well dehned 
physical meaning. The cutoff gives the range of the validity of the model. There are 
a few implementations in four dimensional theories: sharp momentum cutoff in 3 or 4 
dimensions, modihed operator regularization (based on Schwinger proper time method 
[3]). In the Nambu-Jona-Lasinio model diherent regularizations proved to be useful 
calculating diherent physical quantities jl]. 

Using a naive momentum cutoh the symmetries are badly violated. The calculation 
of the QED vacuum polarization function (Il^i,{q)) shows the problems. The Ward 
identity tells us that = 0 , e.g. in 

(i) 

the two coefficients must be the same n(g^). Usually the condition 11(0) = 0 is required 
to dehne a subtraction to keep the photon massless at 1-loop. However this condition 
is ambiguous when one calculates at 7 ^ 0 in QED or in more general models. For 
example in the case of two diherent masses in the loop, it just hxes n(g^, mi, m 2 ) 
in the limit of degenerate masses at = 0 . Ad hoc subtractions does not necessarily 
give satisfactory results. 

There were several proposals how to dehne symmetry preserving cutoh regulariza¬ 
tion. Usual way is to start with a regularization that respects symmetries and hnd 
the connection with momentum cutoh. In case of dimensional regularization already 
Veltman observed [B] that the naive quadratic divergences can be identihed with the 
poles in two dimensions (d=2) besides the usual logarithmic singularities in d=4. This 
idea turned out to be fruitful. Hagiwara et ah [ 6 ] calculated electroweak radiative 
corrections originating from ehective dimension-six operators, and later Harada and 
Yamawaki performed the Wilsonian renormalization group inspired matching of ehec¬ 
tive hadronic held theories |7j. Based on Schwinger’s proper time approach Oleszczuk 
proposed the operator regularization method j8], and showed that it can be formu¬ 
lated as a smooth momentum cutoh respecting gauge symmetries 019]. A momentum 
cutoh is dehned in the proper time approach in PI with the identihcation under loop 
integrals 

kfj,hi/ }■ ~q^uk ( 2 ) 

insteac^of the standard d = 4. The degree of the divergence determines d in the result; 
A^ goes with d = 2 and In(A^) with d = 4. This way the authors get correctly the 
divergent parts, they checked them in the QED vacuum polarization function and in 
the phenomenological chiral model. 

^In what follows we denote the metric tensor by both in Minkowski and Euclidean space. 
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Various authors formulated consistency conditions to maintain gauge invariance 
during the evaluation of divergent loop integrals. Finite m or inhnite iniDS] num¬ 
ber of new regulator terms added to the propagators a’la Pauli-Villars, the integrals 
are tamed to have at most logarithmic singularities and become tractable. Pauli- 
Villars regularization technique were applied with subtractions to gauge invariant and 
chiral models [Til [Til Uni ITT] . Differential renormalization can be modihed to fulhll 
consistency conditions automatically, it is called constrained differential renormaliza¬ 
tion [18]. Another method, later proved to be equivalent with the previous one |19] . 
is called implicit regularization, a recursive identity (similar to Taylor expansion) is 
applied and all the dependence on the external momentum (q) is moved to hnite inte¬ 
grals. The divergent integrals contain only the loop momentum, thus universal local 
counter terms can cancel the potentially dangerous symmetry violating contributions 
[20112T] . A strictly four dimensional approach to quantum held theory is proposed in 
[22] . They interpreted ultraviolet divergencies as a natural separation between physical 
and non-physical degrees of freedom providing gauge invariant and cutoff independent 
loop integrals, it was also applied to non-renormalizable theories [23]. Gauge invari¬ 
ant regularization is implemented in exact renormalization group method providing 
a cutoff without gauge hxing in [23]. Introducing a multiplicative regulator in the 
d-dimensional integral, the integrals are calculable in the original dimension with the 
tools of DREG [25] . 

In this chapter we give a dehnite method in four dimensions to use a well de- 
hned momentum cutoff. We show that there is a tension between naive application of 
Lorentz symmetry and gauge invariance. The core of the problem is that contraction 
with cannot necessarily be interchanged with the integration in divergent cases. 
The proper handling of the terms in divergent loop-integrals solves the problems 
of momentum cutoff regularizations. Working in strictly four dimensions we use the 
conditions of respecting symmetries to dehne the integrals with free Lorentz indices. 
Using our method loop calculations can be reduced to scalar integrals and those can be 
evaluated with a sharp momentum cutoff. We give a simple and well dehned algorithm 
to have unambiguous hnite and inhnite terms [26] dubbed as improved momentum cut- 
oh regularization. The method was successfully applied earlier to a non-renormalizable 
theory [271128] . The results respect gauge (chiral and other) symmetries and the hnite 
terms agree with the result of DREG. There were various other proposals to mod¬ 
ify the calculation with momentum cutoh to respect Lorentz and gauge symmetries 

0 El HH El ED- 

An ideal application of the improved cutoh is the unambiguous calculation of the 
triangle anomaly in four dimensions presented in [29]. Dimensional regularization [1] 
respects Lorentz and gauge symmetries, but as it modihes the number of dimensions (at 
least in the loops) it is not directly applicable to chiral theories, such as the standard 
model or to supersymmetric theories. Gontinuation of 75 to dimensions d 7 ^ 4 goes 
with a 75 not anticommuting with the extra elements of gamma matrices, and it leads 
to “spurious anomalies”, see [30113111321133], and references therein. The loop integrals 
using the novel improved momentum cutoh regularization are invariant to the shift 
of the loop momentum, therefore the usual derivation of the ABJ triangle anomaly 


3 


would fail in this case. We extend the method to graphs involving 75 , and show that 
the proper handling of the trace of 75 and six gamma matrices provides the correct 
anomaly, the {75, 7 /^} anticommutator does not vanish in special cases under divergent 
loop integrals. 

The rest is organized as follows. In section 2 we present how to dehne a momentum 
cutoff using the method of DREG, then we give the gauge symmetry preserving condi¬ 
tions emerging during the calculation of the vacuum polarization amplitude. In section 
4 we discuss the condition of independence of momentum routing in loop diagrams. 
Section 5 shows that gauge invariance and freedom of shift in the loop momentum 
have the same root. Next we show that the conditions are related to vanishing surface 
terms. In section 7 we give a dehnition of the new regularization method and in section 
8 as an example we present the calculation of a general vacuum polarization function 
at 1-loop. In section 9 we show that the QED Ward-Takahashi identity holds at hnite 
order using the new method. Section 10 deals with the famous triangle anomaly and 
the chapter is closed with conclusions. 


2 Momentum cutoff via dimensional regularization 


DREG is very efficient and popular, because it preserves gauge and Lorentz symme¬ 
tries. Performing standard steps the integrals are evaluated in d = 4 — 2e dimension. 
Generally the loop momentum integral is Wick rotated and with a Feynman parameter 
(x) the denominators are combined, then the order of x and momentum integrals are 
changed. Shifting the loop momentum does not generate surface terms and it leads to 
spherically symmetric denominator, terms linear in the momentum are dropped and 
(j 2 |) is used. Singularities are identihed as 1/e poles, naive power counting shows that 
these are the logarithmic divergences of the theory. In DREG quadratic or higher 
divergences are set identically to zero. However Veltman noticed [5] that quadratic 
divergences can be calculated in d = 2 — 2(e — 1) in the limit e —)■ 1. This observation 
led to a cutoff regularization based on DREG. 

Garefully calculating the one and two point Passarino-Veltman functions in DREG 
and in 4-momentum cutoff the divergences can be matched as PE] 


47rju^ ^ + 1^ - 

(3) 

■jE + In (47r + 1 = liiA^, 

(4) 


where /r is the mass-scale of dimensional regularization and je is the Euler-Macheroni 
constant appearing always together with 1/e. The hnite part of a divergent quantity 
is dehned as 


/finite lim 




/(£) - m 


0 


7 e + Indvr 1 ] - R(l) 


e-1 


+ 1 


(5) 


where i?(0), R(l) are the residues of the poles at e = 0, 1 respectively. Note that 
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Fig. 1. 1-loop vacuum polarization diagram 


in the usual e —)■ 0 limit the left hand side (LHS) of ([2D vanishes and no quadratic 
divergence appears in the original DREG. 

The identihcations above dehne a momentum cutoff calculation based on the sym¬ 
metry preserving DREG formulae. This cutoff regularization is well dehned, but still 
relies on DREG. Let us see the main properties in the calculation of the vacuum po¬ 
larization function. In 11^,^ the quadratic divergence is partly coming from a term 
via 2 ■ Qktvk'^i which is evaluated at d = 2 instead of the d = 4 in the naive cutoff 
calculation. The terms cancel if and only if this term is evaluated at d = 2. This is 
a warning that the usual 

k^ku ■^gfiuk'^ ( 6 ) 

substitution during the naive cutoff calculation of divergent integrals might be too 
naive, especially as an intermediate step, the Wick rotation is legal only for hnite 
integrals. A further hnite term additional to the logarithmic singularity is coming 
from the well known expansion in — i (7 + |); cind it is essential to retain 

gauge invariance. We stress that the shift of the loop momentum is allowed in DREG, 
an improved cutoff regularization should inherit it. In the next sections we derive 
consistency conditions for general regularizations. 


3 Consistency conditions - gauge invariance 


Galculation in a gauge theory ought to preserve gauge symmetries. Gonsider the QED 
vacuum polarization function with massive electrons. We start generally (see Fig. 1.) 
with two fermions with different masses in the loop m and restrict it to QED later. 




-Hgf 


d'^k / /c + mg /j + m^ \ 

{2tiY ^ k'^ — {k + qY — nil) ' 


(7) 


is calculated with the standard technique, only the k^ky terms are considered 
with care. After performing the trace, Wick rotating and introducing the Feynman 
x-parameter the loop momentum is shifted {kE^ + xqE^j 


= g^ 



dx 


d^E ^IeYeu - g^u Ye + 7^) - 2x(l - x)qE^qEy + 2x(l - x)g^yq% 
(2vr)4 (/| + A)^ 

( 8 ) 
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where A = x{l — x)q‘^ + (1 — x)rril + xml- In QED rria = mi, = m and 5 ^ = e, it 
simplifies to Ai = x(l — x)q^ + Having a symmetric denominator and symmetric 
volume of integration the terms linear in Ie^l are dropped. After changing the order of 
momentum- and x-integration the loop momentum is shifted with x-dependent values, 
xqe/i and sum up the results during the integration. Different shifts sums up to a 
meaningful result only if the shift does not modify the value of the momentum integral 
(it will be discussed in the next section). 

In QED the Ward identity tells us, that 


= 0. (9) 

In (jH]) the terms proportional to qE fulfill the Ward-identity ([9]) and what remains is 
the condition of gauge invariance 


r', ( d'lE IeM _ 1 /' . /■ dHE 1 

A ‘'D (2ir)y(| + A.f + 


( 10 ) 


This condition appeared already in [HI US\. Any gauge invariant regulator should 
fulfill flTU]) . It holds in dimensional regularization and in the momentum cutoff based 
on DREG of Section 2 . In [TTl 113] a similar relation defined the finite or infinite 
Pauli-Villars terms to maintain gauge invariance. 

So far the x integrals were not performed. Expanding the denominator in the 
x-integration can be easily done and we arrive at a condition for gauge invariance at 
each order of g^. At order g^” we get (omitting the factor (27r)^) 


dH, 


^EfJ'Eu 


(/| + 2n 


" 9^^-^ I d If 


(ll + m^y 


n = 1 , 2 ,... 


( 11 ) 


The conditions (ITT]) are valid for arbitrary mass, so it holds for any function A 
independent of the loop momentum in 1 -loop two or n-point functions with arbitrary 
masses in the propagators. These conditions mean that in any gauge invariant regular¬ 
ization the two sides of flTT]) should give the same result. We will use this condition to 
define the LHS of (ITTl) in the new improved cutoff regularization. This is the novelty 
of our regularizations method. 


4 Consistency conditions - momentum routing 

Evaluating any loops in QFT one encounters the problem of momentum routing. The 
choice of the internal momenta should not affect the result of the loop calculation. The 
simplest example is the 2-point function. In (jT]) there is a loop momentum /c, and the 
external momentum g (see Fig. 1.) is put on one line {k + q, k), but any partition of the 
external momentum {k + q + p, k + p) must be as good as the original. The arbitrary 
shift of the loop momentum should not change the physics. This independence of the 
choice of the internal momentum gives a conditions. We will impose it on a very simple 
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loop integral 


d^k 


J^2 _ ^2 


— / d'^k- 


k^ + p^j, 


= 0 


( 12 ) 


{k + p) — m? 

which turns up during the calculation of the 2-point function. Expanding flT^ in 
powers of p we get a series of condition, meaningful at p, p® .... At linear order we 
arrive at 


d^k 


Pm 


^2 _ ^2 


k^k-p 
2 — m?)^ 


= 0 , 


(13) 


which is equivalent to flTTD for n = 1. At order p^ a linear combination of two conditions 
should vanish 


PpPaP/3 d^k 


4:kak/3 


9aP 


\{k‘^ — m‘^)^ {k"^ — rn?y 


9pp - 


‘Ikakf^kp 9a0kp 

^ ^ (k? — m?y {k‘^ — m?)^ 


= 0 . 


(14) 

These two conditions get separated if the freedom of the shift of the loop momentum 
is considered in J At leading order it provides 

9pu 


P 


d^k 


(fc^ — m?y 


^ kpk„ 

{k‘^ — m?) 


= 0 , 


(15) 


equivalent with flTTl) for n = 2. Using flT5|) twice the second part of the condition (1141) 
connects 4 loop momenta numerators to 2 fc’s. Symmetrizing the indices we get 


^ 4 ^ kakfjkpkp 

{k^-m?f 


1 f ^4:j^9(y-P9pp T 9ap9fip T 9ap9i3p 

24 J [k? — m?y 


(16) 


Invariance of momentum routing provides conditions for symmetry preserving regu¬ 
larization and these conditions are equivalent with the conditions coming from gauge 


invariance. 


5 Gauge invariance and loop momentum shift 

We show at one loop level that gauge invariance of the vacuum polarization function is 
equivalent to invariance of a special loop integrand against shifting the loop momentum 
m- Consider dehned in ([7]), performing the trace we get 


= -g 


2 , kp {k^ + q^) + k^ {kp + qp) - 9 pu ik"^ + k ■ q - mamb) 


(27r)"^ {k‘^ — my) {{k + qY — ml) 

Specially in QED = mb = m, gauge invariance requires ([2]), which simplihes to 

d^k ( kp^qp kp 


(17) 


iq'^Up^iq) = g 


(27r)^ \{{k + q)^ — mY {k"^ — m^ 


= 0 . 


(18) 


This example shows that the Ward identity is fulhlled only if the shift of the loop 
momentum does not change the value of the integral, like in (IT^ . 
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In [21] based on the general diagrammatic proof of gauge invariance it is shown 
that the Ward identity is fulhlled if the difference of a general n-point loop and its 
shifted version vanishes 



We interpret f[T8|) and flT^ as a necessary condition for gauge invariant regularizations. 

6 Consistency conditions - vanishing surface terms 

All the previous conditions are related to the volume integral of a total derivative 



( 20 ) 


The total derivative on the LHS leads to surface terms ra, which vanish for hnite 
valued integrals and should vanish for symmetry preserving regularization. In our 
improved regularization this will follow from new dehnitions. The LHS is in connection 
with an inhnitesimal shift of the loop momentum k, it should be zero if the integral 
of the term in the delimiter is invariant against the shift of the loop momentum. The 
vanishing of this surface terms reproduces on the RHS the previous conditions (IT^ 
and fITT]) . In (120]) starting with any odd number of fc’s in the numerator we end up 
with some conditions, three fc’s for n = 3 provide ([TB]) after some algebra. Starting 
with even number of fc^’s in the numerator on the LHS in 0201 ) we get relations between 
odd number of /c^’s in the numerators, which vanish separately. 

These surface terms all vanish in DREG and give the basis of DREG respecting 
Lorentz and gauge symmetries. Vanishing of the surface term is inherited to any 
regularization, like improved momentum cutoff, if the identihcation f[T0|) is understood 
to evaluate integrals involving even number of free Lorentz indices, e.g. numerators 
alike k^k^. The value of integrals with odd number of /c’s in the numerator are similarly 
dictated by symmetry, these are required to vanish by the symmetry of the integration 
volume. 

7 Improved momentum cutoff regularization 

We propose a new symmetry preserving regularization based on 4-dimensional momen¬ 
tum cutoff. During this improved momentum cutoff regularization method a simple 
sharp momentum cutoff is introduced to calculate the divergent scalar integrals in the 
end. The evaluation of loop-integrals starts with the usual Wick rotation, Feynman 
parametrization and loop momentum shift. The only crucial modihcation is that the 
potentially symmetry violating loop integrals containing explicitly the loop momenta 
with free Lorentz indices are calculated with the identihcation 
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1 


dHi 




m + A)"+' ^ 2n®'"'i ’‘‘‘•^(1% + A)" 

under the loop integrals or with more momenta using the condition ([TB]) or generaliza¬ 
tions of it, like 


74 ? ^Efi^Ei/^Ep^Ea 


-)■ 


QpuQpcr + 9pp9 ua “ 1 “ 9p(t9i'P 

4n(n — 1) 


dH, 


{1% + A) 


n—1 * 


( 22 ) 


The momentum integrals containing further the loop momentum with indices summed 
up (e.g. l\) in the numerator are simplihed in a standard way cancelling a factor in 
the denominator 


dH ® 


IfJ 


Ep^Eu 


’ {ll + A) 


n+l 



• • • 

(‘I + A)" 



((| + A)”+‘' 


(23) 


Integrals with odd number of the loop momenta vanish identically. These identihca- 
tions guarantee gauge invariance and freedom of shift in the loop momentum. Under 
any regularized momentum integrals the identihcations fl2T|) or generalizations like fl22|) 
are understood as a part of the regularization procedure for n = 1, 2,... For hnite inte¬ 
grals (non divergent, for high enough n) the standard calculation automatically fulhlls 
fl21ll22p . The connection with the standard substitution of free indices is discussed in 
Appendix A. 

Fulhlling the condition (ITTli via the substitution fl?Ili the results of momentum cutoff 
based on DREG of section 2 are completely reproduced performing the calculation in 
the physical dimensions d = 4 PZIE5]. The next three examples show that the new 
regularization provides a robust framework for calculating loop integrals and respects 
symmetries. 


8 Vacuum polarization function 

As an example let us calculate the vacuum polarization function of Fig. 1. in a general 
gauge theory with fermion masses rria, mf,. Performing the calculation in 4 dimensions 
generally the Ward identities (required by the theory) are restored by ambiguous and 
ad hoc subtractions. The hnite terms of different calculations do not match each other 
in the literature, see [36], papers citing it and EH. For sake of simplicity we consider 
only vector couplings. Performing the trace in ([7j) we get ([T7|) . Now we can introduce 
a Feynman x-parameter, shift the loop momentum and get ([8]) after dropping the 
linear terms. Generally we are interested in low energy observables like the precision 
electroweak parameters and need the hrst few terms in the power series of 11^,^ (g). 
Using the rule fl2T|) for n = 1 and expanding the denominator in the scalar loop and 
x-integrals can be easily calculated with a 4-dimensional momentum cutoff (A). The 
result in this construction is automatically transverse 

Ilpuiq) = £3 {q^9py - qpQu) [n(0) + g^n'(O) + ...] . (24) 
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The terms independent of the cutoff completely agree with the results of DREG j35] 
the logarithmic singularity can be matched with the 1/e terms using (jl]). Up to O 
we get 

n(0) = jiml + ml) mb)^ In f - (25) 

4 2 \mambj 

m\ +ml - 2mamb {ml + ml) \ 

4 (m2-m2) 

The first derivative is 


n'(o) 


2 4m2m2 — Smamb {ml + ml) 1 ^ / A2 \ 

9 6(m2-m^)^ 3 ^ \mamb J 

{ml + ml) {ml — imlmf + ml) + Gmlmf ^ f ^l\ 
<o{ml-mlf \^l) 


(26) 


The photon remains massless in QED, as in the limit, m^ = mb we get n(0) = 0. 

The proposed regularization is robust and gives the same result if the calculation 
is organized in a different way. Introducing Feynman parameters and shifting the loop 
momentum can be avoided if we need only the first few terms in the Taylor expansion 
of q. For small q the second denominator in f[T7|) can be Taylor expanded, for simplicity 
we give the expanded integrand for equal masses, up to O(g^) 


^^,u{q) = -9 


d^kf 


2k ^ki/ 


4 


(2^)4 (^1 + ^2)2 

2 {kEiiQEu + kEuQEbi) ■ Qe 


(fcl + m 2 ) 
1 


+ 


4 {kE ■ qEY 

{kY + m^Y 


(27) 


{k^Y + m 2 )'' 


S'Ail' 


4 


2 {kE ■ qEY 

{k‘Y + m^Y {kl; + m 2 )'^ + m 2 )" 


+ 


Taking into account that kE-qE = kEadEa, and fl22]) can be used and the remaining 
scalar integrals can be easily calculated. The result agrees with (12^ and (126]) and the 
finite terms with DREG if and only if we use the proposed symmetry preserving sub¬ 
stitutions. Applying the naive kE^ikEu —>■ ^gfiukE substitution ([6]) in both approaches 
the finite terms will differ from each other and also from the result of DREG. This is 
why finite terms differ from each other in [2^ and |37) . 

The gauge invariance of the improved momentum cutoff regularization can be 
checked directly by the well known identity in QED. 


9 The Ward-Takahashi identity 

In this section we show by explicit calculation that the QED Ward-Takahashi identity 
is fulhlled for infinite and finite terms using the proposed regularization at 1-loop. The 
identity reflects the gauge invariance of the underlying theory. 
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Fig. 2. 1-loop diagrams for the Ward identity 


Following the notation of |38) it has to be proved that 

= -5T>^ip,p) 


dp^ 


'^=m 


(28) 


j/=m 


where S is the electron self-energy (see Fig. 2. left panel) 


iu{p)Tiu{p) = —e^ / dx 


dH _ —2x -|- 4m 

7 — ^uip) - oUip) , 

(27r)‘ '^^’(P-A^ + iff 


(29) 


here A 2 = —x{l — x)p^ -|- (1 — x)m^ + xp?, I = k — xp, m is the mass of the electron 
and /i is the infrared regulator. 


dS 


dpp 


y=m 


a 


dx 


-xj^ ( In 


A 2 


1 + 


2 ( 2 -x)(l-a:) 


(1 — x^m? + xp? J (1 — xYm? + xp^ 

dF^ is the electron vertex correction (see Fig. 2. right panel) 

-( bdF^ r i = 2 - 2 / u{p') [/(7^/^+ A) + - 2 m(fc + {k + q)Y] u{p) 

U[p)0 U[p) y + 

After using the Dirac equation in the limit p = p' and g = 0 we get 


(30) 
. (31) 


iu{p)5T^u{p) = 2e^ / dxdydz5{x + y + z — 1) 


dH 


X 


u{j)) / + (x^ — 4x -t- l)m^ 7 '^] u{p) 

(p - A3 + tff ^ 


(32) 


where A 3 = (1 — z^)rn? + zp‘^ and I = k — zp. Here / = 21^1'^'yy — 7 ^/^, for the 
hrst term (|2T1) should be used for n = 2 or directly fl 6 Tl) from Appendix B. After the 
momentum and x, y integration 


a 


^ 


dz 


(1 — 2 ;) ( In 


A 2 


(1 — zYm? + zp^ 


-1 + 


(1 — 4x - 1 - z^) 
(1 — z)m? + zp^ 


(33) 
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Fig. 3. Feynman graphs contributing to the triangle anomaly, k 2 = k, ki = k — qi, 
ks = k + q 2 j ki = k - q 2 and k^ = k + qi. 


The result of the new method is the constant —1 after the log, with the naive calculation 
using fl50p one would get —1/2. Calculating the Feynman-parameter integral taking 
care of the infrared regulator the identity fl 2 ^ holds up to /h? terms at 1 -loop 


dpf. 




%i=m 


a 


j/=m 


/A 2 \ 


= — -In — +ln +2 


27r \2 \m?) 






* A2 ’ A2 


(34) 

We have seen that the proposed method provides regularized 1-loop electron self-energy 
and vertex correction in QED which fulhll the Ward-Takakashi identity. 


10 Triangle anomaly 


In the improved momentum cutoff framework the triangle anomaly has to be recalcu¬ 
lated. The loop integrals using the novel improved momentum cutoff regularization 
are invariant to the shift of the loop momentum, therefore the usual derivation of the 
ABJ triangle anomaly would fail in this case (cannot pick up a finite term shifting the 
linear divergence). We extend our regularization to graphs involving 75 , and show that 
the proper handling of the trace of 75 and six gamma matrices provides the correct 
anomaly, the { 75 , 7 ^} anticommutator does not vanish under divergent loop integrals. 
In this section we show that the new method provides a well dehned result for the 
famous triangle anomaly. 

Consider the 1-loop triangle graph on the left on Fig. 3. 


rjnfllPp _ 


^ ( 5 + m ^J^ + m ^ /c+ /^2 + m \ 

(27r)‘^ (^k - qP'^ - m'^'^ k"^ - m?'^ [k - rrp'^ } ' 


(35) 


The amplitude of the crossed graph is similar with (gi, /i) and (^ 2 , J^) interchanged 
(rp^ivp _ rppvp Ward identities require 



0 , 

(36) 

= 

0 , 

(37) 

(gi + g2)pT^"" = 


(38) 
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where corresponds to the same graphs with a pseudoscalar current instead of the 
axialvector one. There is a formal proof of fl35]) . Replace 

- (gi + q2)pYl^ = - (/:+ - m) 7® + (/;- ^ - m) 7^ ( 39 ) 

The first term combines with the numerator of the last term in fl35|) and cancels the 
denominator. If 

{ 7 '‘, 7=}=0 ( 40 ) 

assumed, then the second term in ( 159 ]) is — {/k— /h — m) 7® = +7^ (/:— /h — m) + 
2m7^. Here the first term cancels the adjacent fraction in ( 155 ]) and the second term 
gives the right hand side of ( 158 ]) . The {—{qi + q 2 )pTi’^^ — difference is 

2 f d^k ( 5 /ii + m ^Jt + m ^ 5 ^Ji + m ^ /i+ ^2 + m \ 

"" J (27r)4 y (fc - ^ k^- (k + q^f -mV' 

(41) 

Shifting k ^ k + qi'm the first term and passing 7^^ through 7^ (using again (l4Up l to the 
back of the second term we arrive to a formula that is totally antisymmetric under the 
interchange of (gi, /i) and (g2, ^^), and thus adding the crossed graph the result 

vanishes. Similarly (156]) and (157]) can be proven but in this case (I50]) is not needed to 
apply, because the terms leading to cancellation are not separated by a factor of 7^. 
The loop momentum can be shifted, this is a fundamental property of the improved 
momentum cutoff regularization. 

However (I361I38P cannot be all true. Pauli-Villars regularization or careful simple 
momentum cutoff calculation identifies a finite anomaly term when shifting the linearly 
divergent integral. There is still a remaining ambiguity in connection with momentum 
routing and which identity contains the anomaly term in (I361l38p . At the same time 
in iimroved momentum cutoff or DREG (1361) and (]57p holds but the proof of (13 8 p is 
falsqj, it relies additionally on (ITUp . This is the first sign that the naive anticommutator 
(irni) cannot be used in all situations. 

The explicit calculation of the triangle diagram (155p is based on the evaluation 
of the trace of 7® with six 7’s. There are various methods to calculate this trace 
with superficially different terms at the end. The different results of the trace can be 
transformed to each other using the Schouten identity, involving two loop momenta it 
reads 

k (^pu\p T k kfid-ai/Xp T kjjk dpaXp T kxk dpuap T kpk 6pi/Xa 0. (42) 

In the present method this identity cannot be used for the loop momentum [k) of a 
divergent integral before applying the identifications (l2Tp or (|22P , because it would mix 
free Lorentz indices and contracted indices, which must be evaluated in a different way 
(DREG faces the same difficulty). After performing the identifications (12ip and (I55p 
the quadratic loop momenta factors cancel with the denominators. The remaining 
formula contains the loop momentum in the numerators at maximum linearly, the 
corresponding Schouten identity can be applied. The root of the problem is that in 

^Functional integral derivation of the anomaly shows that the Ward identity corresponding to the 
axial vector current (I55P must be anomalous |39| . 
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case of divergent integrals the totally antisymmetric tensor e^^xp cannot taken out of 
the integral, similarly to the case of g^i, in the previous section. No such problem 
emerges for hnite integrals. 

The breakdown of the early application of the Schouten identity forces us to choose 
one dedicated calculation of the trace. The trace is calculated not using the anticom¬ 
mutator only 

{ip, Iv} = ‘^9pv , (43) 

and general properties of the trace. The unambiguous result is 

— Tr p'l>\ ^ap/ 3 u 9 p\ (^ap 0 p 9 u\ T ^apup 90 \ (^a 0 vp 9 p\+ 

~^^p0up9a\ (^Xap09pv T ^\apu9P0 ^\a0u9pp T ^\p0i'9poi ^\pap9i'0 T 

~^^\pap9vp ^Xpp09uQ T ^Xpua9p0 ^Xpup9a0 T ^Xpu09oip- (44) 

It reflects the complete Lorentz structure of the 7 matrices in the trace. This choice 
of the trace also appeared in earlier papers without detailed argumentations gnisi]. 
All different calculations of the trace are in agreement with each other and with if 
fHnj) is modihed. 75 and 7 ^ does not always anticommute (rather the anticommutator 
picks up terms proportional to the sum of a few Schouten identities.) Explicitly, the 
following dehnition will eliminate all the ambiguities burdening the calculation of the 
trace of 75 and six 7 ’s 

Tr [{7p, 75 } ixlalplplu] = 2Tr [^^p757A7a7p7/3 - ^/3p757A7a7p7v+ 

+5'pp757A7a7/37^^ - 9apl5'lx'lplplv + 5'Ap757a7p7/37^^] • (45) 

The above anticommutator is defined only under the trace. 0451) can be understood as 
the { 75 , 7 p} anticommutator is dehned by picking up all the terms when moving 7 p all 
the way round through the other hve 7 ’s. Evaluating the trace the right hand side is 
proportional to Schouten identities. Under a divergent loop integral it will not vanish in 
the present method (nor in DREG). The nontrivial anticommutator contributes to the 
triangle anomaly but vanishes in non-divergent cases and for less 7 ’s. The amplitude 
of the triangle diagrams can be calculated with the dehnition of the trace fl44p and the 
identihcations fl 2 T]) . fl 2 ^ . Finally we arrive at the extra anomaly term in fl38l) . 

In what follows we calculate directly the anomaly term missing in fl38j) . We use 
fl3^ and move (/:— /ii — m) from the back to the front in fl35l) using flT3l) . Without 
this trick the same result is obtained evaluating the trace of six 7 ’s and 7 ^ using fl44)) . 
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which is consistent with non-anticomniuting 7 ^ in this special case, see 


“ (?i + 0.2)pTi'^^ — 


d^k 

1 


Tr 


-7 


/ii — m jk — m 


-|- 7^7^-7 


/c + m /:+ /i 2 + 'm 


+ 


+ 27 ® 


- 27 ® 


+ 27 ^ 


- 27 ^ 


, (fc-gi)(fc + g 2 ) 

/k— /ii — m ' /k — m' {k + ^ 2 )^ — 


-7 


7 


1 


1 


- ql) 

/k— /ii — m ' Jk — m ' Jc+ /i 2 + rn 


-7 


r 


+ 




1 


{k - qi){k) 


/k— Jii — m Ji+ J [2 + m k"^ — m? 

1 1 ^ {k^^-ql) 

-^ - 

jk— — m Jc + m Jc+ /12 + m 


(46) 


With algebraic manipulations using the antisymmetry of the trace including 75 and 
four 7 ’s we can group the terms 


“ (?i + q2)pTi'^^ — 


d^k 

W? 


Tr 7 ® 


Jk , J^ A2YY , o 2 / 4 i A2YY , 

H---T-h 2m^————-h 


N 1 N 2 


NoN 


2i'l3 


N 1 N 2 N 3 


- Ai A2l^Y ■k‘^+ A ^2^'!'' ■ kqi + 


WiV2A^3 

+ A ■ kq 2 + A A AY k^+ A A2Y Ak"'} + 

+ .r .r ( + A ^lY Y ' ^2 “ A A7^7'' ' 


WA^2A^3 

- A A2YY ■ A- A A A 7 Af- A AiY A A 


(47) 


where Ni = {{k — qiY — m^), W = {Y — m?) and W = {{k + ^ 2 )^ — rn?)- The first 
two terms vanish after performing the trace and the integral (they are proportional 
qiaqiji and e^^“^g 2 aQ' 2 « respectively). The third one gives 2m times the pseu¬ 
doscalar amplitude = TY'' + , 




r dY 

1 

1 

_N,N2N3_ 


(48) 


we get interchanging (gi,/i) -f-)- (^ 2 , in the integrand. 

The last five terms in fITTj) contain one factor of the loop momentum {k) and 
after tracing vanish by the Schouten identity, the loop integration does not spoil the 
cancellation. The contribution of the one but last five terms in the curly bracket does 
not vanish. It contains two factor of the loop momentum, and it is proportional to the 
Schouten identity fl42li broken under the divergent loop integral. Calculating it with 
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the improved momentum cutoff of Section 2 using the formulas of the Appendix (or 
with DREG) we get the anomaly term 

- {qi + (49) 

In the case of the naive substitution ([ 6 ]) the Schouten identity fH2l) is satisfied, the 
curly bracket vanishes. (In that case with simple momentum cutoff the anomaly term 
originates from shifting the linearly divergent first two terms in fHTli . but the result 
depends on momentum routing.) The presented method identifies without ambiguity 
the value of the anomaly in the axial-vector current and leaves the vector currents 
anomaly free without any further assumptions. 


11 Conclusions 

We have presented in this chapter a new method for the reliable calculation of di¬ 
vergent 1 -loop diagrams (even involving 75 ) with four dimensional momentum cutoff. 
Various conditions were derived to maintain gauge symmetry, to have the freedom of 
momentum routing or shifting the loop momentum. These conditions were known by 
several authors dH [1211201 El]. Our new proposal is that these conditions will be satis¬ 
fied during the regularization process if terms proportional to loop momenta with even 
number of free Lorentz indices (e.g. ~ k^k^) are calculated according to the special 
identifications fl 2 T]) and fl 22 |) or generalizations thereof. In the end the scalar integrals 
are calculated with a simple momentum cutoff. The calculation is robust - at least at 
1 -loop level - as we have shown via the fermionic contribution to the vacuum polariza¬ 
tion function. The hnite terms agree with the one in dimensional regularization in all 
examples. The connection with DREG is more transparent if one uses alternatively 
the kfj,ki, —)■ or fl55|) substitution and d takes different values determined by the 

degree of divergence in each term fl52l l53l l5T|) . We stress that this new regulariza¬ 
tion stands without DREG as the substitutions fl 2 T]) . fl 22 |) and scalar integration with 
a cutoff are independent of DREG. The success of both regularizations based on the 
property that they fulfill the consistency conditions of gauge invariance and momentum 
shifting. 

At 1-loop the finite terms in the improved momentum cutoff are found to be equiva¬ 
lent with DREG. For practical calculations we propose to use the same renormalization 
scheme, MS or MS subtractions plus BPHZ forest formula as with DREG. DREG is 
not just the generally used method, but it is proved to be a mathematically rigorous 
regularization within the Epstein-Glaser framework m- The equivalence of the results 
of the proposed method and DREG gives a hint that the improved cutoff method with 
e.g. MS subtraction and BPHZ can be used as a renormalization scheme for more 
complicated diagrams. 

Regularization schemes based on consistency conditions have been applied to more 
involved cases. Gonstrained differential renormalization is useful in supersymmetric 
H and non-Abelian gauge theories, it fulfills Slavnov-Taylor identities at one and 
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two loops IS]. Implicit regularization j2Ul |2I] requires the same conditions as we 
used and it was successfully applied to the Nambu-Jona-Lasinio model | 2 U| and to 
higher loop calculations in gauge theory. It was shown that the conditions guarantee 
gauge invariance generally and the Ward identities are fulfilled explicitly in QED at 
two-loop order IZH. In an effective composite Higgs model, the Fermion Condensate 
Model I2S] oblique radiative corrections (S and T parameters) were calculated in 
DREG and with the improved cutoff, too, the finite results completely agree. The 
calculation involved vacuum polarization functions with two different fermion masses 
and no ambiguity appeared |2Z1135]. 

As an application the triangle anomaly was calculated within the 4 dimensional 
improved momentum cutoff framework. The property that the loop-integrals are in¬ 
variant under the shift of the loop momentum spoils the usual derivation of the ABJ 
anomaly in the presence of a cutoff. We calculated the trace fl44p corresponding to the 
triangle graphs of Fig. 3. (75 and six 7 ’s) and the Ward identity (75 and four 7 ’s) 
exploiting only the standard anticommutators of the 7 matrices flT3|) and not using the 
7 ^, 7 ^ anticommuting relation. It turns out that different evaluations of the trace agree 
with each other if and only if 7 ® does not always anticommute with 7 ^, rather { 7 ^, 7 ®} 
picks up terms proportional to the Schouten identity fl45|) if it is multiplied with five 
more 7 ’s under the trace. The trace of the { 7 ^, 7 ^} anticommutator multiplied with 
three 7 ’s vanishes, as the definition of ^( 7 ^ 70 , 7 ^ 7 ^ 71 ,) is unambiguous. The right hand 
side of fl45|) is only non-vanishing if it is under a divergent loop momentum integral. In 
four dimensional field theory the nontrivial properties of 7 ^ and 7 ’s appear first time 
in the divergent triangle diagram. 

Traces involving 7 ^ and even number of 7 ’s can be calculated in the same manner 
avoiding the anticommutation of 7 ^and 7 ^. First the order of 7 ‘"’s are reversed applying 
then using the cyclicity of the trace we get back the original trace in the reversed 
order, the difference gives the trace twice. This way the { 7 ^, 7 ^} anticommutator can 
be defined, it will not vanish generally. If it is multiplied with (2n-l-l) 7 ’s it is equal to 
the sum of (2n+l) trace involving 7 ^ and (2n) 7 ’s, see fl45|) . It is well known that the 
general properties of the trace and { 7 ^, 7 ®} = 0 are in conflict with each other, this led 
to the’t Hooft-Veltman scheme BEQ]. Our proposal similarly modifies { 7 '^, 7 ®} but 
works in four dimensions and the modifications come into action only under divergent 
loop integrals involving enough number of 7 matrices. There were attempts to preserve 
{ 7 ^, 7 ^} = 0, but in that case the cyclicity of the trace was lost jESj. We have shown 
with the new method that the vector currents are conserved and the axial vector current 
is anomalous, and no ambiguity appears. 

The new regularization is advantageous in special loop-calculations where one wants 
to remain in four dimensions, keep the cutoff of the model, like in effective theories, 
in derivation of renormalization group equations, in extra dimensional scenarios or in 
models explicitly depending on the space-time dimensions, like supersymmetric theo¬ 
ries. Similar approaches succeeded in the calculation of the anomalous decay of the 
Higgs boson to two photons in four dimensions, where gauge invariance is crucial 
IISIIIT]. We argue that the method can be successfully used in higher order calcula¬ 
tions containing terms up to quadratic divergences in (non-Abelian) gauge theories. 
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as it allows for shifts in the loop momenta, which guarantees the ’t Hooft identity 
pn 1^ . This symmetry preserving method can be used also in automatized calcula¬ 
tions (similar to [19]) as even the Veltman-Passarino functions j50] can be defined with 
the improved cutoff. The strength of the improved momentum cutoff method is that it 
can be used in theories with quadratic divergencies important for example in gauge the¬ 
ories including gravitational interactions El. The calculation in the Einstein-Maxwell 
system was presented in [52] and quadratic contributions to the photon 2-point func¬ 
tion were identified but after renormalization they vanished and did not change the 
original running of the gauge coupling. 


A Connection with simple momentum cutoff 


What is the relation of the new method with the standard (textbook) 


(50) 


substitution? We have to modify it in case of divergent integrals to respect gauge 
symmetry, i.e to fulhll flTTl) . Lorentz invariance dictates that in flTT]) the LHS must be 
proportional to the only available tensor i.e. 


lE^llEv —t -^gnJE (51) 

can be used, where d is a number to determin^. Now both sides of equation flTT]) can 
be calculated with simple 4-dimensional momentum cutoff. The different powers of A 
can be matched on the two sides, and for n = 1 we get the following conditions (from 
gauge invariance) for the value of d, 




/ A'^ + m‘^\ 

\ y 

1 

d 





— for finite terms. 
4 


(52) 

(53) 

(54) 


We see that for finite valued integrals when the Wick-rotation is legal, the condition 
f[TT|l and the rule fl2T]) gives the usual substitution ([50]), but for divergent cases we get 
back the identification partially found by [610^ and others. Quadratic divergence 
goes with d = 2, logarithmic divergence goes with d = 4 plus a finite term (a shift), it 
is the -|-1 in equation (|1|). For more than 2 even number of indices generalizations of 
fl5ip should be used, for example in case of 4 indices the 


^EfJ'Eu^Ep^Ea t 


d(d + 2) 


i^QpvQpu T 9pp9ua T 9pcr9up) ^E' 


(55) 


^The usual method is to calculate the trace (and get d= 4 ), but interchanging the order of tracing 
(multiplication with g^'^) and calculating the divergent integrals cannot be proven to be valid. 
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substitution works. 

We emphasize again that for non-divergent integrals the rules d^T]) and fl2^ give 
the same result as the usual calculation (150]) . 


B Basic integrals 


In this appendix we list the basic divergent integrals calculated by the regularization 
proposed in this paper. In the following formulae can be any loop momentum (k) 
independent expression depending on the Feynman x parameter, external momenta, 
etc., e.g. A(x, q, rria, rrih). The regulated integrals are denoted by /^^.gg meaning 
the integration is understood for Euclidean momenta with absolute value below A. The 
integrals fl56|) and flbOj) are just given for comparison, those calculated with a simple 
momentum cutoff. 


d^k 


Areg *(27r)^ k"^ — m? (dvr)^ 


In 


A^ -|- m? 




d^k kf^ku 




Areg *( 2 vr)^ (fc 2 

d'^k k^kukpka 


1 9iiv 


A^ — In 


A^ -|- 


Lg ^(2vr)4 {k^ - 
r d'^k kdkpk^ 

Lg ^(2vr)4 (F - 
‘‘ d‘^k 


(4- 


(56) 

(dvr) 2 2 m2 J J ^ ^ 

7r)2 8 V \ w? JJ ^ ’ 

2 A^ — 3m^ In 


1 9pv / r, A 2 


Areg *(2vr)^ (fc2 
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(dvr) 2 4 
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A^ -|- m^ 


m^ 


+ m — 
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A? -|- m? 


■ m- 


2 \ 2 


d^k kpki, 


Areg *(2vr)^ (fc2 

r p4 L 


■ m- 


2\3 


d'^k k^kpk, 


Lg *(2vr)4 (fc2 - ~ (4^12 

r Mu u u u u 1 „ 


(dvr) 2 

1 9ph 
(dvr) 2 4 

1 17/ri' 
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A^ -|- m^ 


m^ 
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A 2 _|_ 


- 1 


In 
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A^ -|- m? 
m 2 

A^ -|- w/ 
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m 


A 2 - 1 - m 2 


+ 5 


m 


m 


d k kpkiikpkfj _ 1 9 pi^9pcr~k9pp9 vo-T 9pcr9vp 


A2 m2 (A2 + m2)" 
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In 


A^-|-m^ 


m^ 


m 


A 2 -|-m 2 


(59) 

(60) 
(61) 

4 ) (62) 

-lj(63) 


fl56tl58p depend on the same function of A. fl571158]) are traced back to fl56|) via 
and fl22|) . fl5^ and fl6^ have a different A dependence. Evaluating these integrals at 
hrst step fl2^ is used, then or fl22|) can be applied to the remaining free indices. 
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